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0\ ■ Abstract 

m _ 

CN ! We present a model independent analysis of new-physics contributions to the rare 

f^^ ! decays K^ -^ tt+z/z/ and Kl ~^ 7r°z/z/. We parameterize the effects of new physics in 

^ ! these decays by two parameters: r^ and the phase 6k, with r^ = 1 and 6*^^: = in 

the Standard Model. We show how these parameters can be extracted from future data 
together with the relevant CKM parameters, in particular the angle (3 of the unitarity 



a^ 



p. 



P-t. 

D i triangle. To this end CP asymmetries in B ^ i'Ks and B -^ tt+tt as well as the ratio 



|Vufe/V^fe| have to be also considered. This analysis offers simultaneously some insight 
in a possible violation of a "golden relation" between K — > vrz/z/ decays and the CP 
rS \ asymmetry in i? ^ i^Ks iii the Standard Model pointed out some time ago. We illustrate 

5^ \ these ideas by considering a general class of supersymmetric models. We find that in 

the "constrained" MSSM, in which 9k = 0, the measurements of 'Bi{K^ -^ tt^uu) and 
Br(i^L -^ Ti^vv) directly determine the angle (3. Moreover, the "golden relation" remains 
unaffected. On the other hand, in general SUSY models with unbroken R-parity the 
present experimental constraints still allow for substantial deviations from rj^ = 1 and 
6k = 0. Typically 0.5 < rj^ < 1.3 and —25° < 6k < 25°. Consequently, in these models 
the violation of the "golden relation" is possible and values for Br(i^'^ -^ -k^vv) and 
Br(i^L -^ TT^uiy) departing from the Standard Model expectations by factors 2-3 cannot 
be excluded. Simultaneously, the extraction of the "true" angle jS from K -^ ttuu is 
not possible without additional information from other decays. Our conclusions differ in 
certain aspects from the ones reached in previous analyses. In particular, we stress the 
possible importance of left-right flavour- violating mass insertions that were not considered 
before. 



1 Introduction 

Among the CP asymmetries in B decays, the CP asymmetry in Bd — ^ "ipKs ^ is unique 
as it is the only one which in the Standard Model measures directly one angle of the 
unitarity triangle, the angle (3, without any hadronic uncertainties. 

On the other hand, it has been pointed out in [^ that a similar role among the K 
decays is played by Kl — ^ n^uu and K~^ — > vr+i/z/. Indeed, taken together, these two 
decays offer a very clean determination of sin 2/3 by measuring their branching ratios only 
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and Pc = 0.40±0.06 represents the internal charm contribution to the amplitude A{K^ —>■ 
TT^ui?) which is known including next-to-leading QCD corrections [^. The error in Pc is 
dominated by the renormalization scale uncertainties and the value of rric. 

Using the well-known expression for sin 2/3 in terms of the time-integrated CP asym- 
metry a^Ks i^ Bd -^ ipKs decay, one finds an interesting connection between rare K 
decays and B physics [|2|: 

2r,(5i,fi2) _ ^ l+xl 
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or equivalently 

[sin 2/3]^^^ = [sin 2/3] ^j^g , (5) 

which must be satisfied in the Standard Model. Here Xd measures the B^ — B^ mixing. 
As stressed in [0, this "golden relation" involves, except for P^? only directly measurable 
quantities. Due to very small theoretical uncertainties in (^), this relation is particularly 
suited for tests of CP violation in the Standard Model and offers a powerful tool to probe 
the physics beyond it. These points have been recently reemphasized in [^ Q. 

The present status of the Standard Model predictions for K — > -kvv can be summarized 
byi: 

Br(K+ -^ TT+z/z/) = (9.1 ± 3.8) - 10"^\ 
Bi{Kl -^ 7i%u) = (2.8 ± 1.7) - 10"i\ (6) 



where the errors are dominated by the uncertainties in the CKM parameters. On the 
experimental side, the first event for K'^ -^ tt^uu has been recently observed by the 
BNL787 collaboration [0, giving 

Br(X+ ^ 7r+z/z7) = (4.2+|^) ■ IQ-^", (7) 

in the ball park of the Standard Model expectations. The most recent upper bound on 
Br(irL -^ vrVz/) from FNAL-E799 is 1.8 ■ 10"^ A new proposal at Brookhaven, AGS2000 
|§], expects to reach the single event sensitivity 2 ■ 10~^^, allowing a 10% measurement of 
Br(i^L —>■ vr°z/z/). With a similar accuracy for Bt{K^ — > tt'^uu) a measurement of sin 2/3 
through (|lD with an error of ±0.05 is possible 0. This is comparable with the expected 
accuracy for sin 2(3 from a^Xs i^^ the first years of the next decade at B factories. 

The purpose of this paper is to investigate how the relation (H) could be affected by 
new physics beyond the Standard Model. In particular we will consider a large class of 
supersymmetric models. This will also allow us to reach some general conclusions on 
supersymmetry effects in Kl —>■ n^ui? and K~^ -^ ti^uu. 

Our paper is organized as follows. In section 2 we briefly recall the derivation of the 
formulae above in the manner useful for generalizations. In section 3 a model independent 
analysis of new-physics effects in i^ ^ ttuu decays is presented. In order to obtain some 
insight in the violation of the relation (^), this analysis is then combined with the model 
independent analysis of the unitarity triangle presented in Q]. In section 4, Kl -^ ir^vv 
and K^ —>■ ir^vv are analyzed in a large class of supersymmetric models. This analysis 
gives some insight in the violation of the relation (H) in these models. There we compare 



our results with a recent analysis of Nir and Worah |T^]. Our conclusions differ in certain 
aspects from the ones reached by these authors. Section 5 summarizes the main results 
of our paper. 



2 The Case of the Standard Model 

It is instructive to recall first how (|l]) is derived. To this end let us consider the effective 
Hamiltonian for K — » t^vv: 



V2^sm ew,Jr^.,^ ^ 



which originates in Z-penguin and box diagrams with internal charm (X^vl) and top quark 
(X(xi)) exchanges. Here Aj = V*.Vid and Xt = m^/M^. The dependence on the charged 
lepton mass resulting from the relevant box-graph is negligible for the top contribution. 
In the charm sector this is the case only for the electron and the muon but not for the 
r-lepton. 



Now Kl -^ Ti^vv is purely CP violating in the Standard Model with CP violation 
proceeding dominantly in the decay amplitude [|ll|]. Let us then introduce the quantity 
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with T]x = 0.985 being the QCD correction calculated in [^. Expressing the matrix 
elements < 7i^\sL'j'^dL\KL > and < 7i^\sL'y'^dL\K'^ > through Br(K+ — > vr'^e"''//) one 
finds: 

and 

Here 

and 

where we have used 
1 



K+ = rK+ - 
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sin^ Gw = 0.23, Bt{K+ -^ 7r°e+z/) = 4.82 ■ 10"^ . 
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Finally, rx+ = 0.901 and r^L = 0.944 summarize isospin breaking corrections in relating 
K^ -^ -K^vv and Ki^ —>■ tt^uu to K^ —>■ ir^e^u respectively [|l^. The known potential 
ambiguity in the value for sin^ Bw can be reduced by including two-loop electroweak con- 
tributions to K ^ T[uu [0. With sin^ Bw = 0.23, these two-loop electroweak corrections 
are estimated to be neghgible: at most 1 — 2%. 



Given 'Bi{K^ -^ tt+z/z/) and Br(ifL -^ ti^vv) one finds: 



ReF = -eiJBi- B, 



ImF = e2\ B2, 



(17) 



with the reduced branching ratios defined in (|^). Si being equal to ± indicate the four-fold 
discrete ambiguity which must be resolved by using other decays. 

In the Standard Model, the present knowledge of the CKM matrix implies ei = + 
and £2 = +• Furthermore X^l and X{xt) are real and the contribution proportional 
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Figure 1: The "true" unitarity triangle. 



to ImAc = — ImAf can be safely neglected in view of X^i^/ X{xt) ~ 0(10 ^). Using the 
Wolfenstein parameterization fl^ we have 



and consequently 

and 

Moreover one finds 



K = -A, Xt = -A^X'^Rte-'P 

ReF = -\Pc + A^RtX{xt) cos^] 
ImF = A^RtX{xt)?,inl3. 

1] = Rt sin P 
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A^x{xty 

where rj is the height of the unitarity triangle and Rt is the length of one of its sides as 
shown in fig. |^. Consequently, the measurements of Br(i^"'" -^ tt^vv) and 'Bi{Ki -^ tt^uu) 
allow the construction of the unitarity triangle. Inserting (|19]) and (^) into (|1^) one 
derives (0) with a = 1. The corrections 0{X^) in (^ follow from the improvement [^ of 
the Wolfenstein parameterization. In order to simplify our presentation of the effects of 
new physics we will neglect these corrections. They can be included in a straightforward 
manner if necessary. 

It should be stressed that ([T2|), (|1^) and (|1^) are rather general as they are valid in 
any extension of the Standard Model, provided 

• the new physics contributions to the tree level decay K^ — > ir^e^u and 

• the contributions of operators different from SL'~f^dLu[^'j^i'[^ 

can be neglected. Consequently, they are useful for generalizations of this discussion to 
models in which F includes additional contributions not present in the Standard Model. 



The neglect of new-physics contributions to K^ —>■ 7r°e+z/ is very well justified in all 
known extensions of the Standard Model. On the other hand, as pointed out in ||^, if 
lepton number is violated also operators Si7^(iLZ/^7^^'|^ may contribute. In this case the 
final state ir^uu is not necessarily a CP eigenstate and CP conserving contributions to 
i^L -^ TT^vv can in principle be substantial. In our opinion such a situation is rather 
unlikely and will not be considered in this paper. 



3 A Model Independent Analysis 

New physics can modify the effective Hamiltonian in (^ through new box diagram and 
penguin diagram contributions involving new particles such as charged Higgs, charginos, 
stops etc. In the case of Kl — > tt'^vv also new contributions to K^ — K^ mixing should in 
principle be considered. However, the smallness of ^i^ = 0(1O~^) implies that the K^ — K^ 
mixing effects should be safely negligible in K^ —>■ -k^vv ||1I| independent of the model 
considered as long as Br(i^L -^ ii^uu) > (9(10^^^). Since the Standard Model prediction 
amounts to Br(i^L — ^ n^uu) = (2.8 ± 1.7) ■ 10^^^ we expect that this inequality is 
satisfied also in a large class of its generalizations. Consequently the only new physics 
contributions relevant for i^L -^ Ti^vv are the ones which modify the Hamiltonian in (^. 
Since this Hamiltonian governs both K^ -^ Tt^vv and Ki^ —>■ tt^uu, the discussion of 
new physics effects is considerably simplified compared to the case in which the K^ — K^ 
mixing had to be considered. This situation is opposite to the one encountered in the CP 
asymmetry a^Ksi where one expects the dominant new contributions from new phases in 
B^ — B^ mixing with negligible contributions from new phases in the decay amplitude. 

In order to simphfy the analysis we will assume that: 

• the unitarity of the three-generation CKM matrix is maintained; 

• the new-physics contributions in semi-leptonic tree-level B decays used to determine 

|\4{,| or A as well as iKifc/Kb] can be neglected. 

Then the new-physics contributions io K ^ txvv can be described quite generally by two 
new parameters r^ and Ok defined by 

X„ew=r,,e-^^^X(xi). (22) 

Here Xnew summarizes all contributions to Ties except for the Standard Model contribu- 
tion with internal charm quarks. In a given model tk and 9k are generally complicated 
functions of masses of new particles and of new couplings. They can also carry additional 
irit dependence resulting from loop diagrams in which simultaneously the top quark and 
new particles are exchanged. Well-known examples of such a situation are box and Z- 
penguin diagrams with top and charged Higgs exchanges present in two Higgs doublet 
models and models based on supersymmetry. 
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In principle we could also modify the first term in (P). This, however, would un- 
necessarily complicate our discussion. Any new-physics effects on the term proportional 
to Ac can be included in the second term without the loss of generality. Yet one has to 
remember that Ac 7^ At and any new physics contribution proportional to Ac will not only 
modify tk but also give a contribution to 6k even if this new contribution does not carry 
any new phase. 

In this context, it should also be remarked that at least in the two Higgs doublet model 
the loop diagrams with charm and charged Higgs exchanges give negligible contribution 
to K'^ -^ -K^vv because of very small charm- if ^ Yukawa couplings. 

With X{xt) replaced by Xnew we have 



and 



which implies 



ReF = -[Pc + A^RtrKX{xt) cos{f3 + Ok)] 
ImF = A^RtrKX{xt) sin(/? + 9k) , 
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P^ + A^RlrlX\xt) + 2PA^RtrKX{xt) cos(/? + 9k)\ (26) 



and 



Br(irL -^ ttVp) = Ki^ ■ A^Rlr\X'^{xt) sm^{l3 + 9k) 
where r^ is given by 



eiV^i -B2-Pc 



£2VB2 

For ei = +, 62 = + and cr = 1, r^ reduces to (H). Finally, 
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A^X{xt) 



RtVK sm{P + 9k) = r]f- 



(27) 
(28) 
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In order to simplify the discussion, we will assume that ei = + and £2 = + as in the 
Standard Model. 

We make the following observations: 

• If the new-physics contributions are governed by the CKM matrix and no contribu- 
tions proportional to Ac as well as no new complex couplings are present {9k = 0) 
then, even with rK 7^ 1, sin 2/? can be directly found by measuring Br(i^+ -^ ti^uu) 
and Bt{Ki^ -^ n^uu) only. It does not require the knowledge of Vk which depends 
on unknown masses of new particles and new couplings. Examples of such a sit- 
uation are the two Higgs doublet model and the simplest supersymmetric models 
such as the "constrained" Minimal Supersymmetric Standard Model (MSSM), to 
be specified in the following section. 
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Figure 2: The "fake" unitarity triangle from K —>■ ttuu. 



Yet even in such situation the numerical value of (3 extracted using eq. (|25|) may 
generally differ from the value one would find using the Standard Model analysis 
of the unitarity triangle, which is based on B^^ ~ B^s mixings and Ek- Indeed, 
new-physics contributions to the latter quantities may require a modified value of 
/3, which as we have seen can be directly extracted from eq. (531) if 6'i^ = 0. 



In general, however, Or ^ ^ and as seen in ( pSf ) Br(if+ -^ n^vu) and Br(i^L ~^ 
Ti^vv) allow to extract only {(3 + 9x) instead of (3 as in the case of the Standard 
Model. Once {13 + Ok) has been determined, one can use (p6D or (|^) to find the 
product RtTR- 



We also observe that the ratio on the left-hand side of 
in (^) but a new quantity ?7f. 



no longer measures rj as 



Thus in the presence of new physics the measurements of 'Bi{K^ -^ tt^vv) and 
Br(i^L —* vr^z/z/) imply a "fake" unitarity triangle shown in fig. ^ in which (3 and 7] 
are replaced by {{3 + 9k) and rji respectively. This implies that the side {Rb)i in this 
fake triangle will generally differ from i?b defined as 



Rh 



K 



ub 



V, 



cb 



(30) 



In order to extract P, 9k, Rt and tk additional input is necessary. 



Let us first remark that in order to see whether 6'^^ 7^ and r^^ 7^ 1 one can simply 
check whether {Rb)i in fig. ^ agrees with R^ in ( |5DD which is expected to be unaffected by 
new-physics contributions. Similarly one can use the standard analysis of the unitarity 
triangle which involves the Standard Model expressions for ek, B^ — B^ mixing and 
iKife/Kbl- This analysis gives the values of /? and Rt, albeit with substantial uncertainties. 
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Inserting these values into (^6]) and (p7|) and comparing with the measured values of 
Bt{K~^ -^ Ti^vv) and Br(i^L -^ n^uu), may give the first hint whether r^ 7^ 1 and 
9k 7^ are indeed required. 

Now the new physics will generally contribute also to Ek and B^ — B^ mixing. Con- 
sequently in order to extract the values of 6^ and r^ a more refined analysis is required. 
To this end one can use the model independent construction of the unitarity triangle 
of Grossman, Nir and Worah |P which involves the CP asymmetries a^xs ^^d a^^r in 
Bd — *• il^Ks and B^ — » tt+tt^ decays respectively. This analysis is based on the following 
plausible assumption: 

• the b — > CCS and b -^ uud decays for a^Ks ci-nd a^^r respectively are mediated by 
Standard Model tree level diagrams with negligible contributions from new physics 
in decay amplitudes. This implies that the only impact of new physics on a^Ks ^-nd 
Otttt will come through B^ — B^ mixing. 

The remaining three assumptions involving K^ — K^ mixing, the dominance of the 
Standard Model contribution in the determination of iKtb/Khl and the unitarity of the 
three generation CKM matrix, made in [§], are the same as we have made in connection 
with K —>■ TTuu. 

A weak point of the analysis of is the use of the asymmetry a-^T^ to extract the 
angle a. It is well known that this extraction is affected by the "QCD penguin pollution". 
The recent CLEO results for penguin dominated decays indicate that this pollution could 
be substantial [|I^. The most popular strategy to deal with this "penguin problem" is 
the isospin analysis of Gronau and London |[T^. It requires however the measurement of 



Br(i?2 —>■ TT^TT^) which is expected to be below 10^^: a very difficult experimental task. 
For this reason several, rather involved, strategies [0 have been proposed which avoid 
the use of B^ -^ ir^n^ in conjunction with the asymmetry a-j^T^. It is to be seen which of 
these methods will eventually allow us to measure a with a respectable precision. In what 
follows we will assume that all these problems will be overcome and a will be measured 
through a,r7r one day. 

Under these circumstances the new-physics contributions in the analysis of [Q can be 
described analogously to (|22|) by two parameters r^ and 6^, defined by 



where Hgg'(Ai? = 2) and T-C^^{AB = 2) denote the full Hamiltonian (including the new 
and Standard Model contributions) and the Standard Model Hamiltonian respectively. 

With new-physics contributions to B^ — B^ mixing the asymmetry a^Ks does no longer 
measure sin 2/5 but sin2(/5 + 6^) and the relation (^) is not satisfied if 6k 7^ 6^. 

[sin 2(/? + dK)Uu + [sin 2(/? + Qa)\^K^ (32) 
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Figure 3: sin2(/3 + 9k) as a function of Bt{K^ 
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3-10 ^^ and 1-10 ^° and two choices of {61,62)- The results for (+, — ) and (— , +) are 
obtained by reversing the sign of sin2(/5 + 6k) in these plots. 

Since 9k originates in new contributions to the decay amplitude K -^ txvv and 9^, in new 
contributions to the B^ — B^ mixing, it is very likely that 9k 7^ 9^- 

Now as demonstrated in ^, the knowledge of a^Ks^ '^■k-k and Ri, allows to extract in a 
model independent manner /?, Rt and 9d and consequently the true unitarity triangle. In 
order to find Td also the B^ — B^ mixing parameter Xd subject to hadronic uncertainties 
has to be considered. Inserting Rt and (3 into (|25|)-(p9D allows then a model independent 
extraction of 9k and tk from K -^ ttuu. This is also evident from fig. ^ This in turn 
gives more insight in the violation of the relation (H) as anticipated in 



Having extracted the values of Vd, 9d, tk and 9k in a model independent manner, one 
can then compare these values with explicit model calculations of these quantities. Such 
a comparison can either exclude certain models or put bounds on their parameters. 

Clearly in a general case one would have to take care of discrete ambiguities repre- 
sented by £1^2 in (|28D and (|29|) . Analogous discrete ambiguities are also present in the 
analysis of ref . . In order to resolve these ambiguities, additional experimental infor- 
mation will be needed. In fig. H we plot sin2(/? + 9k) as a function of Bt{K~^ — > n'^ui?) 



for different values of Br(i^L -^ tt^vv) and the choices (+, +) and 



for {61,62). The 



results for (+, — ) and (— , +) can be obtained by reversing the sign of sin2(/? + 9k) in 
the left and right plot in fig. ^ respectively. The black square is in the ball park of the 
central values expected in the Standard Model. We observe that even a modest accuracy 
for sin2(/? + 9k) should be sufficient to distinguish between various choices for {61,62). 
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4 K ^ Tiuu and Supersymmetry 

We will now study the effective Hamiltonian for s — > dvv transitions in the framework 
of a generalized SUSY extension of the Standard Model with minimal field content and 
unbroken R-parity. 

Supersymmetric contributions to the process K^ — >■ vr+z/z/ were already considered 
in refs. [|^-[H] and, very recently, in ref. [|l^, which also analyzed Kl -^ ir^uu. With 



the exception of refs. fiy, |2^, the phenomenological analyses have been done prior to 
the top quark discovery and using much less stringent constraints on SUSY masses than 
those presently available. This makes the comparison of our numerical results with these 
papers very difficult. Moreover, most of the previous studies were performed in the 
framework of the so-called "constrained" MSSM, in which universality of the soft SUSY 
breaking terms is assumed. In this case, all SUSY contributions to FCNC processes are 
still proportional to the CKM mixing angles. Consequently, due to the heaviness of the 
superpartners, they are generally expected to be small in comparison with the Standard 
Model contributions. Furthermore, in the case of the "constrained" MSSM, the only 
non- negligible SUSY contribution is proportional to Aj, and therefore no new phase 6k 
is introduced in eq. (|^). As already discussed in sec. ||, in this case it is still possible to 
extract sin 2/5 from K -^ vrz/z/ decays, since eq. (|I]) remains unaffected. 

However, it should be noted that the "constrained" MSSM is based on very strong 
assumptions which do not find at present a strong theoretical motivation. It is therefore 
useful to consider a more general definition of the MSSM, in which the assumptions about 
the universality of the soft SUSY breaking terms are relaxed. In this "unconstrained" 
version of the MSSM new sources of flavour and CP violation are present in the mass 
matrices of sfermions, and in general large contributions to FCNC and CP violating 
processes are expected in this case. To be able to study the contributions to s — *> duu 
transitions in such a generalized extension of the Standard Model, and to take properly 
into account the constraints coming from measured low-energy FCNC and CP-violating 
processes, one needs some kind of a model-independent parameterization of the flavour- 
and CP-violating quantities in SUSY. Such a parameterization has been formulated in 
the framework of the mass- insertion approximation |5^], and will be discussed in detail 
below. 



4.1 Computation of SUSY Contributions 

There are several one-loop SUSY contributions to s — > dui? transitions, according to the 
kind of diagram (penguin or box) and to the virtual particles running in the loop: i) 
charged Higgses and up- type quarks (see fig. ^, ii) charginos and up-type squarks (and 
charged sleptons in box diagrams, see fig. ||), iii) gluinos and down-type squarks (see 
fig. H), iv) neutralinos and down type squarks (and sneutrinos in box diagrams, see fig. |^. 
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Figure 4: Feynman diagrams for the charged-Higgs contribution to s — *> dvv transitions. 
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Figure 5: Feynman diagrams for the chargino contribution to s ^ dvv transitions. 
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Figure 6: Feynman diagrams for the gluino contribution to s ^ dvv transitions. 
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Figure 7: Feynman diagrams for the neutrahno contribution to s ^ duu transitions. 
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First of all we note that, due to the left-handed chirality of the neutrinos and due 
to Lorenz invariance, the quark fields entering in dimension-six operators must have the 
same chirality. Therefore, there are only two possible dimension-six local operators for 
s -^ dvv transitions: SL'^^diVi'^^'i'L and SR-ff.dRUL'J^UL. 

Let us consider first gluino- mediated transitions. As already remarked in ref. ||21[| , 



these contributions are negligible. First of all, there are no gluino-mediated box diagrams. 
Concerning the penguin diagrams, it is clear that if the Z coupling to squarks were 
proportional to the photon one, the amplitude would be suppressed at least by a factor 
q/Mz with respect to the Standard Model one, since by gauge invariance there is no 
contribution to the effective sd'-f vertex of zeroth order in the external momenta (we denote 
by q the momentum of the Z boson). In fact, gauge invariance constrains the effective 
sd'y vertex to be of the form FiSi {q^'~i^ — q^</i) di + F2SRa^^q''dL. This is indeed the case 
if the squarks running in the loop have fixed "chirality" . To avoid the q/Mz suppression, 
one has to perform a double chirality flip in the down-type squark propagator, which is 
once again highly suppressed under the plausible assumption that the left-right mixing of 
squarks of the first two generations is negligible^. 

Let us consider next the operators that the contributions in figs. §, ^ and ^ could 
generate in the effective theory. "RL" operators as SRa^yq^di^L'^^i^L are suppressed by 
q/Mz and either by a quark mass or a chirality flip in the down-type squark propagator. 
"RR" operators as s^'-y^d^ui'-^'^i'L can only be generated by a neutralino exchange through 
U(l) gauge couplings and turn out to be negligible. We are left with "LL" operators as 
SLliidi^L'^^T-'L- In a general supersymmetric extension of the Standard Model lepton 
flavour violation can be present due to a misalignment between lepton and slepton mass 
eigenstates. Therefore operator involving different neutrino flavours can be generated. 
However, the experimental limits on /i — i> 67, r — ;> /i7 and other lepton flavour violating 
processes make their contribution to if — >■ nvv negligible. The operators we will consider 
are therefore the same as in the Standard Model. As in that case, we could have a 
dependence of the coefficient on the neutrino flavour involved, here due to a possible non- 
degeneracy of the sleptons running in the boxes. For our purposes, it is sufficient to use a 
"mean" value for the slepton masses of a given charge, without distinguishing among the 
three neutrinos, as done for the top contribution in the Standard Model. 

Our effective Hamiltonian for K -^ -kvv is therefore given by eq. (^ with X{xt) 
replaced by 

X^,^ = X{xt)+XH{xtH) + C^ + CN, (33) 

where XtH = mf/mjj± and the function Xh^Xih) corresponding to the charged Higgs 
contribution is given in the Appendix. C^ and Cn denote respectively the chargino and 
neutralino contributions. 



^This is also the reason why the sin^ Qw part of the Z couphng does not contribute to the total 
amplitude. 
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While the flavour structure of the Standard Model and charged Higgs contributions 
involves just the CKM matrix, Cn and C^ depend on the transition matrices between 
quark and squark mass eigenstates, whose determination requires the diagonalization of 
the 6x6 squark mass matrices. In the so-called super-CKM basis, in which the flavour 
structure of the quark-squark-gaugino vertices is the same as in the quark-quark-gauge 
boson vertices, the squark mass matrices are given by: 

A#2 _ / {^y)LL {mDLR \ M'^ - I ^'^u)lL {^DlR \ , . 

\ {™-d)lb^ \™-d)rr J \ {^u)lb. \™-u)rr J 



where 



(m^)iL = m2^ + (m")V^(3-4sin2ei^)cos2/51 



6 

{rrvl)LR = -;,m"cot/?-^^A" 

{uiDrr = m2^ + (m")' + ^M|sin2ew. cos 2/51 
(ml) LL = m|^ + (m'^) - -^ (s - 2sin2 6^^) cos2/5 1 

/ 9 X rl ^ V COS 3 . . 

(m|))iH = -Atmnan/? Tf"^ 

imDnn = in]^ + [m'^f -^Mlsm^Qwcos2(31. (35) 

We denote by rn|^, ml^; mj and m| the soft-breaking masses for squarks and by A" 
and A*^ the soft-breaking trilinear couplings for up and down squarks in this basis, at the 
electroweak scale. The quark mass matrices are diagonal in this basis: 

m" = diag (m„, nic, nit) , 

m^ = diag {rrid, rris, rrib) . (36) 

In order to obtain simple analytic formulae in a generic model, it is useful not to 
perform the exact diagonalization of the squark mass matrices, but instead use the mass- 
insertion approximation. Since we are interested in s ^ d transitions, we switch to a 
basis in which the d\ — (P^ — N^ and d\ — u'i^ — Xn couplings are flavour diagonal, and the 
flavour change in the left-handed sector is exhibited by the non-diagonality of the sfermion 
propagators. Denoting by A the off-diagonal terms in the sfermion mass matrices, the 
sfermion propagators can be expanded as a series in terms of 5 = A/m^, where fh is the 
average sfermion mass. As long as A is significantly smaller than m^, we can just take the 
first term of this expansion and compute any given process in terms of these 5's0. This 

^Care must be taken in the truncation of the expansion at the first order in those models where there 
is a strong hierarchy between different mass insertions. For example, in the constrained MSSM, in the 
notations of eq. (||), {ml,)d^si^/m'^ ^ ('7i|,)d,,fc,.(TO|,)6^si/m'*. 
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is equivalent to a first order diagonalization of the squark mass matrices around their 
diagonal part. 

The basis we have chosen is obtained starting from the super-CKM basis by rotating 
the left-handed up-type squark fields: 



ul 



V^ul, (37) 



where V is the CKM matrix. In our basis, the up-squark mass matrix is given in terms 
of eqs. (|3^) and (^) by 

with the notation 

(i^DundL = {^DrlV {m1f)unun^ = {mDnR. (39) 

In the above formulae, u and d represent the three-generation space: u stands for u, c, t 
and analogously for d. 

In principle one should diagonalize exactly the up-squark mass matrix in the ti — Ir 
sector before performing the mass-insertion expansion. However, this would make the 
equations more cumbersome without adding any new feature in the results, and therefore 
we will stick to the basis above. It is a straightforward exercise to include this exact 
diagonalization in our formulae. 

The chargino (C^) and neutralino {Cn) contributions are given by 

Cn = XiyRg^dLi (40) 

where the functions X^ and X^, which depend on SUSY masses and respectively on 
chargino and neutralino mixing angles, are given in the Appendix. The R parameters are 
defined in terms of the mass insertions in the following way: 

r>D _ K^^DJsLdL tjU _ \'^u)sLdL 

-Tie, rl, ^ ,.^ , r, -fl. 



where m| and m|^ are respectively the squared masses of the down- and up-type squarks 
of the first two generations. 

Eq. (^OD shows explicitly the dependence of the various contributions on the off- 
diagonal entries of the squark mass matrices, that are completely unknown apart from 
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quantity 









R 



u 

'LtR 






upper limit 



Min <^ 231 



37 



'-112-55i) 



-112 -54i) 



500 GeV 



500 GeV 



sooGevy ,43^xe^^0<</><2vr 



(5^fev)'xe^^O<0<2vr 



Table 1: Upper limits for the R parameters. Notice that the phase of R^^ir ^^'^ ^^ cLl ^^ 
unconstrained. See Section [4.2.2| for details on how these upper bounds are derived. 



upper limits on their real and imaginary parts, and are strongly model dependent. The 
dependence on masses and other supersymmetric parameters is on the other hand con- 
tained in the X functions, whose order of magnitude is rather model independent and 
can be read from a scatter plot for random values of the relevant SUSY parameters (see 
below) . 

The model dependence is therefore in this way almost completely contained in the R 
parameters. We have chosen the normalization of the Rs in such a manner that their 
absolute value in the constrained MSSM is expected to be of order one or smaller. Indeed, 
a rough estimate gives 

3 



(-Rsidi)MSSM 



and 



{R^ 



^t^jMSSM 



(-Rsidi)MSSM 



{RtRdL)MSSM '^ 



Ml 



^y^og 



(47r)2 
At + ^ cot (3 



m. 



(42) 



(43) 



UL 



where Yt is the top Yukawa coupling, YtAt is the (3,3) element of the trilinear scalar cou- 
phng in the up sector and Mq is the universality scale. For Mq = O(Mgut) 



(-^SLdL)MSSM 



Oil) and therefore (i?f^dL)MSSM 

rameters are discussed below and summarized in table ^ 



. log i^Y? -- 
1. Experimental limits on the R pa- 



(47r): 



4.2 Phenomenological analysis 



We are now ready to write a formula for r/^e 



-iB,, 



rxe 



-i0K 



X 



H 



X' 



1 H H + 

XX 



yLL 



yLR 
^X jdU I 



^X pU 






D 

Shdh 



(44) 
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quantity 


range 


tan/3 


2^6 


fi 


(-300 ^ 300) GeV 


M, 


(100 ^ 300) GeV 


^in 


(150 ^ 300) GeV 


"^Ql 


(200 ^ 500) GeV 


^L, 


(100 ^ 300) GeV 



Table 2: Ranges for the SUSY parameters used in the scatter plots. Too hght charginos 
are rejected. 



where X = X{xt). The amount of SUSY corrections to K -^ ttuu branching ratios and 
sin 2/3 depends on how much tkg''^^'^ deviates from unity. This depends of course on the 
size of the mass insertions, that is largely unknown. 

4.2.1 Results in the MSSM 

In figures ||, ^ |iy and |11] the ratios Xl^/X and Xj^ /X are reported as scatter plots for 
random values of SUSY parameters in the ranges given in table |^. The ratio XujX is 
given separately in fig. ^ as a function of mu for tan/5 = 2 (it scales as l/tan^/5). 

Whereas figs. | and 12 are general in the sense that their contributions to r^e"*^^ 
do not depend on the values of mass insertions, the other ones have to be multiplied by 
the corresponding R parameter to get the contribution to r^^e"*^^. On the other hand, 
an upper limit for the contribution of the "constrained" MSSM, and of any other SUSY 
extension of the SM in which all the R parameters are 0(1), can be directly read from 
these figures. Therefore in this quite large class of models the neutralino contributions are 
negligible. The only possible non-negligible contributions come from the charged Higgs 
exchange, in the "corner" of parameter space in which the Higgs mass is near its lower 
limit and tan/5 is very close to one, and from the term X^ in chargino exchange. It is 
interesting to note that these two contributions have opposite signs and therefore tend to 
cancel each other, making the effects very small. This cancellation is particularly effective 
if one takes into account the constraints coming from h -^ s'~f, which imply a correlation 
between the mass of charginos and charged Higgses. We are therefore able to conclude 
that in this class of SUSY models in which the R parameters are 0(1), the effects on r^ 
are certainly smaller than 10%. Moreover, if the R parameters are real, as it is the case 
in the "constrained" MSSM, we have 6k = 0. 

At this point, we would like to comment on the smallness of the chargino and neu- 
tralino contributions, assuming the absence of anomalously large mass insertions. In the 
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Figure 8: The ratio X^jX as a function of \\i\. 

limit of large squark masses, the box diagrams are suppressed by M^jmh relative to pen- 
guin diagrams because the 1/M| coming from the Z propagator in the penguin diagrams 
is replaced by a 1/m? coming from the additional superparticle propagator. One could 
then think that the penguin diagrams might still give contributions not suppressed in 
that limit. Actually this is not the case, once again because of gauge invariance, which 
requires the presence of at least one gaugino-higgsino flip in the chargino/neutralino line 
and an overall Mf/m? suppression as in the box-diagram case. 

4.2.2 Results in a general SUSY extension of the Standard Model 

We now turn to the discussion of the effects which can be obtained in a general SUSY 
extension of the SM. In this case, we do not assume any particular model, but instead let 
the R parameters vary within the range allowed by the available experimental constraints. 

Our numerical analysis uses the constraints on off-diagonal mass insertions found 
in ref. ||25|. To obtain these constraints, it has been imposed that the gluino- mediated 
contributions to K^ — K^ ^ D^ — Z)° and B^ — B^ mixing, as well as to the decay h -^ s'j, 
proportional to each single off-diagonal mass insertion, do not exceed the experimental 
value. Barring accidental cancellations, this analysis gives the largest possible value for 
the R parameters and consequently for the quantities as 9k or r^ considered in this paper. 
On the other hand, in order to simplify our analysis, we have kept the CKM parameters 
at their central values (see table |). In principle, also these parameters should be varied 
in the allowed ranges and this variation should be correlated with the variation of the 
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Figure 9: The ratio \X^^\/X as a function of |/i| 
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Figure 10: The ratio \X^ \/X as a function of |/i|. 
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Figure 11: The ratio \X]\f\/X as a function of |/x|. 
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Figure 12: The ratio Xh/X as a function of niH (GeV), for tan/? = 2 (it scales as 
l/tan^/3). 
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quantity 


value 


Vtd 

Vts 


0.0076- 0.0029 i 
-0.039- 0.0006 i 



Table 3: Values of the relevant CKM matrix elements used in the numerical analysis. 



R parameters in such a manner that the experimental data on all the quantities listed 
above are reproduced. Such a very involved numerical analysis is beyond the scope of this 
paper. Consequently, the numerical results presented for this general SUSY extension of 
the Standard Model should be considered only as order-of-magnitude estimates. 

We now discuss in more detail the constraints on the off-diagonal mass insertions 
which are relevant in our case. From eqs. ( ^OD and (0), we see that the mass insertions 
involved in our expressions are (m|,)^^d^, (m^)^^^^, (m^)^^^^ and (m^)^^^^. 

Let us start with {rn'jj)s^dL. From ref. [^ we immediately get, for degenerate squarks 
and gluinos. 



Re^^^^^^ < 0.04 
rrir 

mi 

dh 



m 



dL 



500 GeV ' 
500 GeV 



(45) 



which implies the upper limit on /2f ^ reported in table |I|. 

We now turn to the constraint on (m^)^^^^. First of all, following ref. |2^, we note 



that SU(2) invariance of the soft-breaking scalar mass matrices implies the following 
relation: 



mi 



V^ 



^ky- 



(46) 



On the other hand, we see from eqs. (|38|) and ( p9D that 



[m 



UJSLdL 



V\ 



m 



UJLL 



V 



21 



rrij + V^' (m 



D) 



V 



V^ 



21 



m^i + (m 



./ + 



6 



(?> - Asm^ Qw) cos2[nV 



where we have used eqs. 
the following constraint: 

Re 



and 



>u) 



si4t 



^II 



<0.04 



. Therefore, we see that, using eq. 
ReV*m',Vtd 



ImilMfi^ < 0.003 



m. 



Ul 



500 GeV 
500 GeV 



+ 



+ 






lmV*m'fVtd 



m. 



Ul 



J 21 

(47) 
we can set 



(45 
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This point was overlooked in ref. [0, in which the authors did not use eq. (|46| ) to set a 
constraint on ("^^)sid£• Instead, they expanded the product V"^ {mjj) llV and used the 
available constraints on (m^)^,^. This yields a weaker constraint on {jn\j)sj^iiLi roughly a 
factor of four larger than the one we obtain. We will return to this point below. 

Concerning the LR and RL mass insertions, things are much more involved. In this 
case, we cannot use the SU(2) symmetry to relate the mass terms in the up sector to the 
ones in the down sector, and the limits in the up sector obtained studying gluino-mediated 



processes are available only for 



ra 



UJcAUBy 



with A, B = L, R. On the other hand, we need 



to know {rnfj)tj^dL and {jnlj)s^tii- These mass matrix elements could be constrained by 
studying chargino contributions to K^ — K^ mixing, B^ — B^ mixing and 6 ^ S7 in the 
mass-insertion approximation. Unfortunately, there is no such analysis available in the 
literature, except for 6 ^ 57 [^, and it would go beyond the scope of this paper to 
discuss these constraints in detail. However, an order-of-magnitude estimate of the limits 
on these matrix elements can be easily obtained in the following way. 

Let us consider K^ — K^ mixing. The gluino contribution to the coefficient of the 

dhl^SL dhl^iSL operator can be written as 



Ml ^ 



(m 



l2 






(49) 



where Cg is a dimensionless function of the average squark and gluino masses, while the 
chargino contribution includes the following terms: 



A. 



{Os,tnytVtd + Vuyt{nii)t^d, 



(50) 



Here Cy. is again a dimensionless function of the average squark and chargino masses 



and of the chargino mixing matrices. We now assume that a'^/M? Cg 



~ at 



./Ml Cy 



which is appropriate to get an order-of-magnitude estimate of the constraints coming 
from chargino exchange. Comparing the chargino and gluino contributions and using the 
constraint reported in ref. ||25|, we obtain, barring accidental cancellations and interference 
effects. 



[m 



UJSLtR 



YtVu ViYAml 



UJtRdL 



fh? 



fh? 






< 0.04- 



ra 



■SOOGeV ^^^^ 

The above equation shows that, due to the presence of the CKM matrix elements in the 
relevant couplings, essentially no constraint can be derived on the (^^)tjjd^ and {jn'lj)t^sL 
entries from K^ — K^ mixing. Analogously, from Bd — Bd mixing one gets 



{mfj)b^tRYtVtd 



v:Mm'' 



UJtRdL 



fh? 



fh? 



<0.1 



m 



SOOGeV 



(52) 



Therefore, we can set the following approximate limit: 



[f^UJtRdL 



m? 



<0.1 



m 



500 GeV 



(53) 
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Constraints coming from chargino exchange in b 
report the following constraint: 



sj have been studied in ref . ||26[ . They 



m: 



U)tRSL 



m? 



<3 



m 



500 GeVy 



(54) 



Constraints on the relevant left-right mass insertions can be also obtained from re- 
quiring the absence in the potential of charge and colour breaking minima or unbounded 
from below direction^. From ref.p^ we get: 



m 



UJtRSL 



J^UJtRdl^ 



< mt\ 2m? + rrih 



2 ) 



(55) 



where m2 is the coefficient of the |i?^2p term in the superpotential. Taking for example 
m2 — m, one sees that the FCNC constraint on (m^)^^^^ in eq. (|53| ) is always stronger 
than the one in eq. (^) in the range of squark masses we consider. On the other hand, 



the constraint on {rn^jt^iSL ^^ sq. (^) becomes tighter than the FCNC one in eq. (^If ) for 
squark masses larger than about 300 GeV. 

The maximum values for the neutralino and chargino contributions consistent with 
present constraints can be estimated multiplying the scatter plots in figs. §-|ll] by the 
upper limits on the relevant R parameters reported in table [^. 

To provide an order-of-magnitude estimate of the possible size of SUSY contributions 
to tk and 6k in this general extension of the Standard Model, we report in the scatter plots 
in figs. and the sum of chargino and neutralino contributions obtained by varying 
the R parameters between zero and the upper limits listed in table |T|. The remaining 
SUSY parameters are varied in the ranges reported in table ^. 

We observe that typically 



0.5 < tk < 1.3 



and 



- 25° <eK< +25' 



(56) 



although values outside these ranges, even if less probable, cannot be excluded. Inspecting 
eqs. (p5|)-(|2^) we find that departures from the Standard Model expectations for the 
branching ratios in (|^) by factors 2-3 are certainly possible. Furthermore, in view of 
substantial values of Ok-, vanishing or even negative values of sin2(/? + 9k) cannot be 
excluded. We recall that the standard analysis of the unitarity triangle gives roughly 
/3~20°±80 i. 

Finally we would like to compare our results with those of ref. [|l^. In that paper, 
K — > TTz/z/ decays have been considered in various specific supersymmetric models based 
on exact universality, alignment, etc. It has been found that in these special models 
SUSY effects were rather small and the extraction of (3 from K -^ t^vv in the presence of 



^We thank Y. Nir and M. Worah for pointing out to us these constraints, which were overlooked in 
the first version of this work. 
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Figure 13: The value of 6k (in degrees) obtained by varying the R parameters between 
zero and the upper hmits given in table |l|, as a function of |/i| (see the text for details). 



SUSY was possible. On the other hand, the authors of ref. |T^ stressed that in generic 
supersymmetric models large SUSY effects are possible. We agree with this statement, yet 
we disagree with the manner this conclusion has been reached. Nir and Worah attribute 
these possibly large effects to flavour violation coming from left-left mass insertions. As we 
have demonstrated in (|48|), a constraint coming from gluino exchange in i^° — K^ mixing. 



not considered in [|10], excludes such large effects from these insertions. On the other 
hand we find that flavour violation in left-right mass insertions, which was not considered 
in ref. ||10[, can give large contributions to K ^ ttuu, even when available constraints are 
taken into account. Another question is related to the contribution of box diagrams with 



chargino exchanges which have not been considered in ref. ||T0[ . In our opinion the neglect 
of these contributions to X^^ cannot be justified. 



4.3 SUSY effects in B ^ ^K^ decays 



We have seen in Sec. T?l that, while in the "constrained" MSSM we still expect to be able 
to measure sin 2/3 from K —>■ vrz/z/ decays with no additional uncertainty with respect to 
the Standard Model case, in a general SUSY model large contributions to 9k are possible 
and therefore the extraction of the true angle /3 from K -^ ttuu becomes impossible. Now 
we would like to briefly comment on SUSY effects in the measurement of sin 2j3 from CP 
asymmetries in B ^ ipKg decays. To this end we use the parameterization (|3TD. 

Let us consider first the case of the "constrained" MSSM. In this model, no new 
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Figure 14: The value of r^ obtained by varying the R parameters between zero and the 
upper hmits given in table |l], as a function of |/i| (see the text for details). 

phase is present {9d = 0) and therefore the CP asymmetry a^px still measures sin 2/?. 
Consequently, in the "constrained" MSSM the relations (|) and (^ are expected to hold. 
However, as we already explained in Section |^, the true value of sin 2/5 that one can obtain 
from a^Ks ^ind K — > irui? will generally differ from the one obtained using the Standard 
Model analysis of the unitarity triangle, which is based on B^^^ — B^^ mixings and Ek 

If we consider instead a general SUSY model, then large contributions to B'^ — B^ 
mixing and, as we have seen, to i^ -^ ttz/z/ decays are possible [^, ^. Moreover, the 
flavour- changing mass insertions that enter in these processes are to a large extent uncor- 
related. This means that large violations of the relations (^ and (^) are possible in this 
case. 



5 Summary and Conclusions 

In this paper we have presented a model independent analysis of new-physics contributions 
to the very clean rare decays K^ -^ tt^uu and Ki^ -^ n^uu. We have illustrated this 
analysis by considering a large class of supersymmetric models. In particular, we have 
investigated whether the relation (^) between K — > vrz/z/ decays and the CP asymmetry 
a^Ks proposed in could be violated in these models. 

The model independent analysis oi K -^ ttz/z/ decays can be formulated in terms of 
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two parameters r^ and 6k in analogy to the parameters r^ and 6^ introduced previously 
in the literature in connection with B^^ — 5° mixing. We have demonstrated in Section 
^ how the parameters r^ and 6k can be extracted from future data on Bt{K^ -^ vr+z/z/) 
and Br(i^L ~^ ii^vv). To this end the model independent analysis of the unitarity triangle 
of ref. P has to be simultaneously invoked. 

Analyzing K -^ vrz/z/ in a large class of supersymmetric models by means of the 



mass-insertion approximation of ref. pj] we arrive at the following conclusions: 



In the "constrained" MSSM, 6d = 6k = 0, implying that the relation @ is well 
satisfied in this supersymmetric scenario. However, the value of sin 2/? extracted 
from K -^ -nvv and a^Ks might differ from the one obtained using the Standard 
Model analysis of the unitarity triangle. 

On the other hand, we have demonstrated explicitly that in a more general class of 
supersymmetric models large deviations from the Standard Model values 6k = ^ and 
rx = 1 are possible. Typically 0.5 < tk < 1-3 and —25° < 6k < +25° implying that 
substantial violations of the relation @) and a departure from the Standard Model 
expectations for 'Bi{K^ -^ tt^uj?) and Br[Ki^ —>■ rr^uu) by factors 2-3 are certainly 
possible. Simultaneously, the extraction of the CKM parameters, in particular of 
the angle /3, from K -^ ttuu alone is no longer possible and a more complicated 
analysis involving other decays is necessary, as described in Section 0. We point 



out that, although we agree with ref. [IC] that in general supersymmetric extensions 



of the Standard Model the effects on tk and 6k can be large, we disagree on the 



calculation of these effects. In fact, the authors of ref. |TD[ state that flavour violation 



in the mass matrices of left-handed up-type squarks could generate contributions to 
K —>■ TTz/z/ as large as the Standard Model ones. However, as we already noticed in 
Section [4.2.2|, there is a constraint coming from gluino exchange in K^ — K^ mixing 



which was not considered in ref. [0; taking this constraint into account, the left- 



left contribution turns out to be small, about 10% of the Standard Model one. 
Furthermore, the contribution to 6k is even smaller. This can be immediately seen 
multiplying the results in figure ^ by the upper limit for -R^^^ reported in table 0. 
On the other hand, we find that flavour violation in the left-right mass matrices for 



up-type squarks, which was neglected in ref. |]TU[, can give contributions to -ft' ^ -kvv 
that are much larger, even when the available constraints are taken into account. 
Thus, large SUSY effects in tk and especially in 6k are mostly due to left-right 
flavour- changing mass insertions, not to left-left ones. 

Independently of what kind of new physics is realized in nature, it is clear that the 
theoretically clean K — >■ vrz/z/ decays in conjunction with the very clean CP asymmetry 
a.^Ks have great potentiality to discover new physics and to show us what this new physics 
is. In particular, a substantial violation of the golden relation (HI) would not only require 
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a generalization of the Standard Model but would also tell us that it cannot be the 
"constrained" MSSM. 
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Appendix 

We collect here the expressions for various functions which appeared in Sect. ^. The for- 
mulas in terms of mass insertions that are shown below are obtained from the general ones 
by using the following observation. Consider a n x n Hermitian matrix A = A^ + A^ with 
A^ = diag(a5, . . . , a°). If the unitary matrix U diagonalizes Ahj A = Wdiag{ai, . . . , an)U 
and / is an arbitrary function, at the first order in A^ we have 

UlfiakW,, = 5.,/(a°) + 4F(a°, aj), (57) 

where 

p/ ^ fix) - f{y) .^„^ 

F{x,y) = . (58) 

x-y 

If the matrix A has degenerate eigenvalues, the limit for x ^ y must be performed 
in eq. ( |55D and one obtains the derivative of f{x). In our case A is the up or down 
squark mass matrix, A^ is its diagonal part and / is the appropriate loop function. We 
are performing an expansion at the lowest order in A^ . Note that in this way no 0{1) 
uncertainties associated with the use of a "mean mass" for squarks of different generations 
are present. Therefore, we can take into account the effects of a light stop. 

We consider the first two generations to be almost degenerate in each squark sector. 
Each term in eqs. (|33D and (^) has a contribution from penguin and box diagrams, except 
for Xh where the box diagram contribution can be safely neglected. 
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We are using the following expressions for the chargino and neutralino mass matrices: 

_ / M2 v^Mvi/sin/3\ 

■ ~ V \/2Mwcosl3 fi ) ' 

(Ml -M^sinOvKCOs/^ M^sinGvi/ sin/3 \ 

M2 M^cos6vi/cos/9 — M^ cos9p|/ sin/3 
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\ M^sinGn^ sin/3 — M^- cos6iy sin/3 
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/ 



In the chargino mass matrix rows 1 and 2 correspond respectively to the W and h^ 
current eigenstates, and columns 1 and 2 correspond respectively to the W^ and h^ 
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current eigenstates. In the neutralino mass matrix rows and columns 1, 2, 3 and 4 
correspond respectively to the B, W^, /i° and /i° current eigenstates. 

The mass eigenstates are obtained by diagonalizing the mass matrices, using 

M^ = K^M^H , Mn = U^M^U, (61) 

with H, K, U unitary matrices and M^ = diag(M^^)„=i^2, M^^ > and M^ = 
diag(M;vJ„=i,...,4, M^„ > 0. 

In the neutralino sector we used the notation 

Unix) = t3{x)UnW3 + taU Qwyix)UnB, (62) 

where x is a generic particle, t^lx) its third isospin component, y{x) is its hypercharge, 
and U is the neutralino mixing matrix defined in eq. (|61|). 

Finally, the loop functions are: 

.. ^ j{x)-jiy) 

j{x,y) = 

x-y 
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k[x, y, z) = k{x, y, z, t) = . (63) 

x-y x-y 
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